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Abstract 



In the framework of loop quantum cosmology anomaly free quantizations of 
the Hamiltonian constraint for Bianchi class A, locally rotationally symmetric and 
J> ' isotropic models are given. Basic ideas of the construction in (non-symmetric) loop 

^ ■ quantum gravity can be used, but there are also further inputs because the special 

Q ! structure of symmetric models has to be respected by operators. In particular, the 

QO basic building blocks of the homogeneous models are point holonomies rather than 

holonomies necessitating a new regularization procedure. In this respect, our con- 
, struction is applicable also for other (non-homogeneous) symmetric models, e.g. the 

' spherically symmetric one. 

u 

^ 1 Introduction 

^ ! Loop quantum gravity a program for the canonical quantization of general relativity, 
^ I has lead to a well-understood theory of quantum geometry predicting discrete spectra 
for geometrical operators like area or volume 0, |], |^, |]. The main open problems 
concern the semiclassical limit and the dynamics of the quantum theory. Recently, a 
strategy to attack the first problem has been proposed by constructing coherent states 
which are peaked around a given classical configuration. 

As for dynamics, there are a number of anomaly free quantizations of the Hamiltonian 
constraint [|l^, [ill and a class of very special solutions . It is, however, not clear how to 



interpret any possible solution, in part because of conceptual problems, most importantly 
the problem of time, which are already present in the classical canonical formulation. These 
classical conceptual questions have been mainly investigated in reduced models subject to 
a symmetry condition such that one has a good control over all solutions. Most prominent 
in this context are cosmological models which are homogeneous in space and present exam- 
ples of mini-superspaces which after reduction can be quantized similar to a conventional 
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mechanical system |TB|, |T1|. In the first part of this series a way was proposed how to 
perform such a symmetry reduction at the kinematical level of loop quantum gravity based 
on the general concept of symmetric states in diffeomorphism invariant quantum theories of 
connections |T^. Equivalently, this can be seen as a loop quantization of a mini-superspace, 
but the quantum states can always be interpreted as generalized, symmetric states of loop 
quantum gravity. Based on the kinematical level, the next step is to quantize the reduced 
Hamiltonian constraint of a given model resulting in the Wheeler-DeWitt operator which 
governs its dynamics. Having such a quantization one can then look whether the classical 
conceptions regarding the problem of time apply and what they imply for the quantum 



dynamics |18 



Quantizations of the Hamiltonian constraints for some models are the subject of the 
present paper. Because the aim is to understand the dynamics of loop quantum gravity, 
we will follow the basic construction in Ref. [|10] as close as possible. Although the en- 
suing Wheeler-DeWitt operators will be quite similar to the one of the full theory, there 
is additional input needed in their quantization, essentially due to the following basic in- 
gredient: in order to obtain the curvature components present in the classical constraints, 
one has to lay a loop in the space manifold whose holonomy approximates the curvature 
of the Ashtekar connection. In the full theory, one then couples the classical limit with 
a continuum limit in which all these infinitesimal loops exactly reproduce the curvature 
components. Moreover, one has to use diffeomorphism invariance in order to ensure that 
the continuum limit does not affect the quantum constraint. In homogeneous models we 
then have to face the following problems: First, there is no place for such loops because 
the reduced models are formulated in a single point and their quantum states are based 
on point holonomies rather than ordinary holonomies. Therefore, a new regularization is 
needed which also provides a substitute for the continuum limit which played an impor- 
tant role in the quantization of the Hamiltonian constraint in the full theory. We will do 
this by extending the point, in which the reduced model is formulated, to an auxiliary 
manifold using the structure of the symmetry group. This step has already been done in 



the construction of derivative operators ||T5[. The next problem then is to lay the loops 



in such a way that the Wheeler-DeWitt operator respects the symmetry, resulting in the 
main difference to the operator in the full theory [ p!0|] . 

We will start by reviewing the structure of loop quantum cosmological models in Section 
and the construction of the Wheeler-DeWitt operator in the full theory in Section |^. Then 
we will be ready to quantize the Hamiltonian constraint of Bianchi class A models (Section 
^ and subsequently of locally rotationally symmetric and isotropic models (Section]^). 

2 Bianchi class A, LRS and Isotropic Models 

In this section we recall the results of a symmetry reduction for cosmological models and 
their quantization at the kinematical level. For details we refer to Ref. . 

Bianchi models are homogeneous models which have a symmetry group S with structure 
constants cfj acting freely and transitively on the space manifold S. Bianchi class A 
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models are characterized by Cjj = 0. Some of them can be reduced by requiring a non- 
trivial isotropy group which is isomorphic either to U{1) (locally rotationally symmetric, 
LRS models) or to SO{3) (isotropic models). In all cases, is a semidirect product of a 
translation group and an isotropy group F. We are interested in the reduction of general 
relativity formulated in real Ashtekar variables and, therefore, need the general expressions 
of invariant SU (2)-connections and CSU (2)-valued density weighted dreibein fields. 

An invariant connection can always be decomposed as A = (p^urWi using the generators 
Tj — —^CTj (with Pauli matrices aj) oi G — SU{2) and left invariant one-forms uj^ on the 
translation group A^ (which is identical to S for Bianchi models). Invariant connections 
are parameterized by the components of a linear map 0: CN CG which is arbitrary 
for general Bianchi models, but restricted to be of the form 

0i = 2-i(aAi + 6A'j , 4 = 2-i(-6Ai + aA'j , cA's (1) 

for LRS models and 0} = cA} for isotropic models. A} is a dreibein which is rotated by 
the internal SU (2)-gauge transformations. 

Similarly, an invariant density weighted dreibein field can be decomposed as = 
^/goplXf with left invariant vector fields on A^ satisfying oo^Xj) = Sj and the determinant 
Qq of a left invariant metric on A^ defined by A u"^ A = ^/god^x. Again, for Bianchi 
models the coefficients p[ are arbitrary, but restricted to be of the form 

pl^2-^paA\+PbAi) , pf^2-^-pkA{+PaAi) , P-^PcAi 

for LRS models and pi = pA^ for isotropic models. 
The symplectic structure is given by 

Wi,p'^} = Ka]5i 

for Bianchi models, 

{a,pa} = {h.,Pb} = {c,Pc} = 1^1^' 

and vanishing in all other cases for LRS models, and {c, p} = kl' for isotropic models. 
Here, k = SttG is the gravitational constant and t' tV^~^ the reduced Immirzi parameter 
with Vo := Jj, d^xy^. 

Reduced expressions can be derived by inserting invariant connections and dreibein 
fields into the unreduced expressions of the full theory. In this paper we need the Euclidean 
parts of the reduced Hamiltonian constraints, which are 

% ri = Qo eijkFjjEjEf, = -eijkCjjCpKPjPk + ^ijkeiim(Pi<PjPjPk 

= -^vkcfAW.pi + ^'A]pi - Mp]pi (2) 

for Bianchi models, 

^-i^(E) ^ _^nm^n{2)^^apa + hp,)p,-n^^^c{pl + pl) 

+{apa + bpb + cpcf - \{aPa + hpbf - {cpcf + \{apb - bpaf (3) 
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for LRS models, where the constants n*-^^ specify the Bianchi type, and 

= -2(n(^) + + n^^^)cp^ + Q{cpf (4) 

for isotropic models, where n*-^-* = n^^^ = nP''^ = for the isotropic flat model (isotropic 
Bianchi I), and n^^^ = n^'^^ = n^^^ = 1 for the isotropic closed model (isotropic Bianchi IX). 

To check the constraint algebras, we will also need the diffeomorphism constraints, 
which are 

Vj = -{Ki'r'cUKPi (5) 

for Bianchi models, and vanish identically for LRS models (provided that the GauB con- 
straint is solved) and isotropic models. So only for Bianchi models demanding an anomaly 
free representation of the constraint algebra will put restrictions on the quantization of the 
constraints, although not for all types (for Bianchi type I the structure constants vanish, 
and so the diffeomorphism constraint). 

For a non- vanishing diffeomorphism constraint, the constraint algebra can be derived 
from the full theory: the Poisson bracket of a diffeomorphism constraint smeared with a 
shift vector field A^" and a Hamiltonian constraint smeared with lapse function yields 
a Hamiltonian constraint with a lapse function given by the shift vector field and the 
derivatives of the lapse function N . In a reduction to homogeneous configurations, the 
lapse function is a constant, and so its derivatives vanish. Therefore, the diffeomorphism 
constraint and the Hamiltonian constraint Poisson commute. Moreover, as can easily be 
checked by a direct calculation using the Jacobi identity for the structure constants c^, 
the diffeomorphism constraint (|^) commutes with the Euclidean part of the Hamiltonian 
constraint (@). 

We now come to a description of the kinematical sector of loop quantum cosmology 



For Bianchi models, the three scalars 0}rj are exponentiated to point holonomies |T9[ taking 
values in SU{2). The auxiliary Hilbert space Tiaux = L'^{SU(2)^ , dyU^) consists of functions 
of these point holonomies being square integrable with respect to the three-fold copy of Haar 
measure. A basis is given by spin network states which combinatorially can be associated 
to spin networks with three closed edges meeting in a single 6- vertex xq. After extending 
Xq to an auxiliary manifold, they can be identified with spin networks embedded in this 
manifold. The diffeomorphism constraint generates transformations which correspond to 
conjugations in S", e.g. rotations for the Bianchi type IX model. 

The auxiliary Hilbert spaces of LRS and isotropic models can be obtained from that 
of Bianchi models by a further reduction which takes care of the fact that for those mod- 
els point holonomies which are rotated by an element of the isotropy group are gauge 
equivalent: while kinematical quantum states for Bianchi models are given by functions 
/(^i,^2,^3) with hi := exp((/)}rj) G SU(2), for LRS models quantum states are given by 
functions 

/(/ii, exp(f Ag)/^! exp(-f A3), h-^) 

where /13 is a function of A3 := A3rj via /13 = exp(cA3), and for isotropic models by 
functions 

/(exp(fA2)/i3exp(-fA2),exp(-fAi)/i3exp(fAi),/i3) . 
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Therefore, their states can be represented with a reduced number of edges. However, 
the configuration spaces are not subgroups of SU{2Y, but only subsets given by a union 
of conjugacy classes. This implies that quantum states can no longer be represented as 
ordinary spin network states, but as generalized spin network states with insertions (see 
Ref. for a discussion). 



3 The Hamiltonian Constraint in the Full Theory 

For real Ashtekar variables, the Hamiltonian constraint in Lorentzian signature is given by 
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n[N] = I d'xN{det gyn^hE^E'; [F^, - 2(1 + l') E^E^^^KlKi) (6) 

where Fab are the curvature components of the Ashtekar connection and 

Kl:= (det q)--2KabE'' 

are the coefficients of the extrinsic curvature of S. The extrinsic curvature coefficients are 
complicated functions of the phase space variables, but luckily they are contained only 
in the second term of the constraint (M). An important step towards a quantization of 



the Hamiltonian constraint has been done in Refs. [[lOi |23| leading to a first consistent 
regularization of the constraint. An important input was a procedure to treat the second 
term of the constraint containing the curvature coefficients K^^, which has been achieved 
by first writing 

n = 2(det g)-5(l + ,2) tT{[Ka, K^W, E']) - H^''^ (7) 

where 

n^^^ = 2(det tr(F„fe[E^ E^]) (8) 

is the Euclidean constraint. In these formulae the variables are written as >CS'C/(2)-valued, 
e.g. E"" := Efr^. The square brackets and tr denote the commutator and trace in this Lie 
algebra. 

It is then necessary to represent the curvature components Ka = K^Ti as quantum 
operators. The recipe of Ref. starts from the observation that the trace 



K := [ d^x^/d^Kabq''^ = [ d^xK^E^ 
Jt, Jt, 



is the time derivative of the volume (independent of the space-time signature) which can 
in a Hamiltonian formulation be written as Poisson bracket of Hamiltonian and volume: 

K = CtV = -{v,n^^^}. 

When one then uses the important fact that {F^, K} = ||2^ to represent the components 
of the intrinsic curvature as 
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translates all Poisson brackets into commutators, expresses the connection coefficients by 
holonomies and uses a quantization of the volume P, |^ one can trace back the quanti- 
zation of the second term of @ to a quantization of the ffist term. 

Furthermore, in Ref. |]T0[ a quantization of the Euclidean constraint H.^^^ has been 
given which completes the quantization of Ti. Here a key point is that in a diffeomorphism 
invariant context only expressions with the correct density weight can be quantized without 
a background structure. This means that the determinant of the metric must not be 
absorbed into the lapse function in order to have a polynomial constraint because this 
would change the density weight. Instead, the key identity 

AT/ 

(det q)--^ [E", E^Y = e^'^'el = 2e"''"^ = 2(/tO"'e"^^{A^, V} (9) 

is used and quantized by turning the Poisson bracket into a commutator, representing the 
connection coefficients by a holonomy and using the known volume operator V. 

To regularize the constraint, the manifold S is triangulated in a way adapted to the 
graph of a cylindrical function on which the constraint operator is to act. This is done 
in such a manner that for each triple of edges meeting in a vertex v a triangle aij, which 
consists of two pieces Sj and Sj of two edges and a third curve aij connecting their endpoints, 
is chosen as a base of the tetrahedron spanned by Sj, Sj and a piece Sk of the third edge. 
The rest of E is triangulated arbitrarily. Using the above identities and an expansion of 
holonomies for infinitesimal length one can see that the vertex contributions 

summed over all triples of edges meeting in v have the correct classical limit for a local 
contribution of the Euclidean constraint. 

Contained in the classical limit is a continuum limit in which the width of the triangu- 
lation vanishes. In the quantum theory, however, such a continuum limit is trivial provided 
we quantize the constraint on the Hilbert space Tidm where the diffeomorphism constraint 
is solved. In this case two operators which are obtained by a triangulation and a refinement 
of it are identical because their actions on a fixed cylindrical function by construction differ 
only by a diffeomorphism moving the edges aij. 

When acting on a cylindrical function the Euclidean constraint changes the graph due to 
the holonomies to the curves aij in such a way that in the neighborhood of any non-planar 
vertex for each two edges incident there a new edge connecting them is generated. The 
action of the operator is well understood, but there are only some rather trivial explicitly 
known solutions and the complete solution space is not under any control. Furthermore, 
there are a lot of ambiguities in the construction of a particular operator and it is not 
clear which one to use. This question can probably be answered only by investigating 
whether the theory has the correct classical limit. A strategy proposed here is to use 
similar regularizations of the Hamiltonian constraint on appropriate sectors of symmetric 
states for a comparison with the corresponding classical mini-superspaces. 



6 



4 Hamiltonian Constraints for Bianchi Models 



We are now ready to present a quantization of the Hamiltonian constraints for cosmological 
models, starting with Bianchi models. Our construction will follow that of Ref. in 



the full theory, but there are important points where new input is necessary. This is 
the case because, first, our regularization has to respect the symmetry conditions and, 
second, because connection coefficients are coded in point holonomies rather than ordinary 
holonomies. Furthermore, Eq. (0) shows that the constraint contains two different parts: 
a part linear in the scalar fields which via the structure constants depends on the 
Bianchi type, and a second part quadratic in the scalar fields which is independent of the 
type. Therefore, although the original expression does not bear any reference to the 
symmetry type, all the reduced models will have different Wheeler-DeWitt operators and 
different dynamics. 

As already mentioned, we will follow the procedure presented in Ref. ||T0| and recalled 



in Section |]when quantizing the expressions (j^) and on the respective kinematical 
Hilbert spaces. Some steps can immediately be copied, for instance we can restrict our- 
selves to the Euclidean constraints because the additional part in the Lorentzian constraint 
can be treated in complete analogy to the full theory by using the extrinsic curvature K. 
Furthermore, we will also make use of commutators of holonomies with the volume oper- 
ator, i.e. this operator will again play a prominent role. At this place it is fortunate that 
it simplifies in symmetric regimes |^ such that an analysis of the matrix elements of the 



Hamiltonian constraint will be easier. 

What we cannot copy in our reduced models is the use of a triangulation of space and 
the subsequent continuum limit. At first sight, we may seem to be in a better situation 
because we are considering a finite dimensional model and do not have to bother about 
a continuum limit. But in the regularization used in the full theory this limit also serves 
the purpose to express the curvature components F^^ entering in terms of closed loops 
in the following way: As recalled above, one starts with a triangulation of space which is 
adapted to the graph of a spin network the Hamiltonian constraint operator is supposed 
to act on. Forming a triangular loop aij based in v one can compute the holonomy hc^ 
appearing in each vertex contribution discussed in Section ^. In the continuum limit any 
such loop shrinks to v such that the holonomies approximate the curvature components. 
Here one makes use of the expansion 

h{A) = l + \eWFi,n + 0{e^) 

for a holonomy h along an infinitesimal parallelogram spanned by two edges of length e 
with directions and f ° in Euclidean space. This is the final ingredient needed to turn 
all contributions to the classical Hamiltonian constraint into quantum operators. 

The ensuing Wheeler-DeWitt operator then contains the holonomies along all the trian- 
gular loops aij as multiplication operators creating new edges aij. Although the construc- 
tion of these loops depends on the triangulation, the action of the operator is triangulation 
independent once it is considered to act on the kinematical Hilbert space where the diffeo- 
morphism constraint is solved. In this sense, the continuum limit is trivial in the quantum 
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theory because for a finer triangulation the newly created edges are just moved using a 
diffeomorphism thus representing the same diffeomorphism invariant quantum state. 

Obviously, we cannot use these techniques to generate the curvature components in 
a symmetric model: First, the prescription to generate new edges violates our symmetry 
conditions (e.g., in the Bianchi models there are just the three edges obtained by smearing 
the point holonomies, and no edge connecting two of them can be generated). Second, we do 
not have a continuum limit in which course loops would shrink and could be approximated 
using curvature components. 

We will present now alternative techniques which are adapted to the symmetry and 
which make use of the methods developed to deal with point holonomies. As in the 



computation of the derivative operators in Ref. |]T5| the regularizing edges will be important. 
All this will be discussed in detail for the Bianchi models and then used in LRS and isotropic 
models. 



As demonstrated in Ref. [|TOl, we have to provide the Euclidean part (0) of the Hamil- 



tonian constraint with the correct density weight in order to be able to quantize it in 
a background independent manner. We thus divide the earlier expression by ^/det q = 

\J\go\ e'^^^eijKPi PjPkl^ multiply with a constant (due to homogeneity) lapse function N 
and integrate over the manifold S: 

n^^)[N] = - ! d'xN{detq)--^goe''''pip-^Fijk 
Jt. 



with 

This can immediately be integrated to 



n^^)[N] = -VoN{l\e'"^-eLMNpfp^p:!\y'e^"'pip^Fuk 
= -2{Kt')-'Ne''''{<p'k,V}Fjj, 



where we used 



and 



V = Vo\n\e'^''eijKPip'p^ 



{0^, eMNLe'^'pfpfpf} = SKL'e'^'^eMNKpfpf (10) 

which, divided by \/det q is essentially the key identity (P) used in Ref. [l^ to quantize the 
inverted triad components. The expression of Fjj in terms of the scalar fields depends on 
the particular symmetry group, i.e. on the Bianchi type, which will be taken care of when 
we choose the route of loops below. 
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4.1 Auxiliary Manifolds 

But before this we have to discuss our regularization scheme which, as aheady noted, 
has to be different from that in the full theory because we cannot use a triangulation 
as regulator and the continuum limit in connection with diffeomorphism invariance to 
remove it. After integrating over E we arrived above at an expression for the Hamiltonian 
constraint which is completely composed of the reduced field components sitting in the 
single point xq of the reduced manifold B := S/S* = {a;o}. Obviously in this reduced 
description there is no substitute for a triangulation of space which could be used as 
regulator. Note, however, that in order to compute derivative operators in point holonomies 



T5| we already had to smear point holonomies along auxiliary edges. To that end we have. 



in the present context of cosmological models, introduced an auxiliary manifold S/F (a 
suitable compactification of the homogeneous space S/F) in which these edges are to lie, 
where in the construction of both the auxiliary manifold and the edges we made use of the 
structure of S/F as & homogeneous space. This structure can now also be used in order 
to regulate the Hamiltonian constraint, but again we cannot simply copy the procedure of 
the full theory (applied to the 6-vertex lying in the auxiliary manifold) because this would 
spoil the symmetry (we are not allowed to create new edges, but only to retrace existing 
ones) . 

Although we then manage to respect the symmetry, we immediately have to face another 
problem: In order to approximate the curvature component by the attached loops and to 
recover the correct classical limit we have to shrink the loops to the vertex, which in the full 
theory was achieved by the continuum limit which we now do not have at our disposal. Our 
only possibility is to shrink the whole auxiliary manifold and with it the auxiliary edges to 
the point xq. In this limit we will recover the classical limit, whereas the quantum theory 
is independent of the extension of our auxiliary manifold. This is completely analogous 
to the full theory, where the diffeomorphism invariant quantum theory is triangulation 
independent. 

Let us now follow in detail the program outlined above. We have a compact homoge- 
neous auxiliary manifold S/F containing the point xq. Each point of S/F can be reached 
from Xq by following pieces of integral curves to the left invariant vector fields Xj (which 
also have been used to define the auxiliary edges smearing point holonomies). If we always 
normalize the invariant vector fields in such a way that their closed (by construction of 
S/ F) integral curves have a fixed parameter length, then shrinking the auxiliary manifold 
is equivalent to multiplying the vector field by a number e smaller than one. So we can 
describe the shrinking of S/ F to S/ F^ by replacing X/ with eXj where eventually we will 
consider the limit e — > in connection with the classical limit. 

We now illustrate the shrinking procedure for the Bianchi I model on the manifold M^. 
In terms of coordinates x^ adapted to the symmetry (the usual Cartesian coordinates), left 
invariant vector fields are Xj = The manifold is non-compact, but we can com- 
pactify it to a three-torus, which then plays the role of S/ F, by restricting the coordinates 
to < < 1 and identifying the points with x^ = and x^ = 1 for each J. By this 
compactification the integral curves of the vector fields Xj are rendered closed and will 
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be used as auxiliary edges based in a base point xq in the three-torus. Denoting the flow 
generated by a vector field X by $f(X): S/F ^ S/F, we can write 



S/F = {$<3(X3)$t,(X2)*t,(Xi)xo -.heS'^ R/Z} . 



The shrinking manifolds 5"/ F^ are then obtained for fixed ij- intervals by replacing Xi with 
eXi everywhere. Then the three-torus as well as the closed integral curves shrink and, for 
infinitesimal e, holonomies along the integral curves (i.e. the smeared point holonomics) 
are approximated by e times the respective scalar field component {hj = 1 + e0/ + 0(e^)), 
and volume integrals over the auxiliary manifold are approximated by e'^ times the reduced 
quantity defined in the base point. 

4.2 Regularization 

Our previous expression for the Hamiltonian constraint was written down in the reduced 
description on the point Xq. To extend it to the auxiliary manifold we just have to integrate 
it, yielding just a factor of because the constraint is constant on S/F^. We thus have 



as point of departure for the quantization. An integration is mandatory here in order to 
obtain an expression without density weight (which can be defined using scale transforma- 
tions respecting the symmetry). The original, unregulated expression can be recovered as 
(in this expression the limit is trivial which will, however, not be the case if we introduce 
holonomy variables below) 



One factor of e can be used to express e(/)|; as infinitesimal holonomy, whereas the remaining 
square of e will be used to express the curvature components via a holonomy along an 
infinitesimal loop, to which we turn now. 

Because the curvature components, expressed in terms of scalar fields, are the only 
place where the Hamiltonian constraint depends on the symmetry group, we have to make 
use of its structure when expressing the curvature by means of suitably laid loops. This 
can be done most easily by composing the loop of pieces of integral curves to the invariant 
vector fields, i.e. as a geodesic parallelogram in an arbitrary homogeneous metric. E.g., 
in the Bianchi I model we can use a chart with coordinates x^ which is adapted to the 
symmetry in such a way that the (commuting) invariant vector fields are given by Xj = 
Using the flows $((X/) generated by these vector flclds, we can form parallelograms in S/F 
composed of the four edges Ci.t i-> $t(X/)a;o, C2:t ^ $((Xj)ci(l), c^:t ^ <l>_j(X/)c2(l) 
and C4: t H- > $_t(Xj)c3(l), for all 1 </ 7^ J < 3, which is closed because the vector fields 
Xi and Xj commute. We can now use the constant scalar fields 0/ = parameterizing 



Fijk 



(11) 




(12) 
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a homogeneous connection A^^ = (p^p^uj^ to compute the holonomy along the parallelogram: 
exp(0/) exp(0j) exp(-0/) exp(-0j) using, e.g. 

KM) = Vexp f dteM^Ti = Pexp /" dtX^cuf ^^^r^ = exp(0^rO . 

J Cl J Cl 

What remains is to note that point holonomies were defined in just the same way by 
calculating holonomies along auxiliary edges defined as integral curves to an invariant vec- 



tor field [15]. We only have to take care that we run through a complete closed integral 



curve to a single invariant vector field which is possible due to our construction of the 
compactification S/F. There are then two possibilities to shrink the parallelogram: First 
we can fix the homogeneous manifold and only shrink the curves by constraining the pa- 
rameter t above to an interval [0, e]. This is the usual procedure for a regularization in a 
three-dimensional theory, but is not applicable here because we can identify a holonomy 
along an integral curve with a point holonomy only if the curve is closed and based in the 
point Xq. Thus we have to use the second shrinking procedure where both the auxiliary 
manifold and the curves embedded into it shrink by using the vector fields eXj as gener- 
ators of both the auxiliary manifold S/F^: and of their own closed integral curves. Now 
the holonomies are always along closed integral curves and can be identified with point 
holonomies for all values of e, but due to 0(eX/) = €0/ (linearity of the map 0: CN —> CG) 
we now have the e-dependent holonomy 

exp(e0/) exp(e0j) exp(-e07) exp(-e0j) 

= (1 + + + e0j + \^(t)]){l - t4>i + ie20?)(l - e0j + ie^j) + 0{i') 

= l + e2(0,0j-0^0,) + O(e3) 

= 1 + e2e,,fc0i0>'= + 0{e') = 1 + e^Fj^jn + 0{e^) . 

Thus we have a product of point holonomies which, for infinitesimal e, reproduces the 
curvature components for Bianchi I. 

For other Bianchi models with a non-Abelian symmetry group, the last equation does 
not recover the correct curvature components because of the additional term containing the 
structure constants cfj of the symmetry group. But, also due to the non-Abelian nature 
of the symmetry group, the parallelogram constructed above will not be closed, even up to 
the order of because now the vector fields used to define the edges of the parallelogram 
do not commute. As is well-known from differential geometry, the endpoint 04(1) of the 
above parallelogram constructed from the vector fields eXj and eXj is a distance [eXj, eXj] 
away from the origin. So we can close our parallelogram by running through the curve 
C5: t (-^ $_(([eX/, eXj])c4(l) after following the curves ci to C4. The ensuing curve is closed 
to the order of e'^ and can be closed by connecting the points 05(1) and xq which affects 
the holonomy along the loop only up to the same order. As holonomy we obtain 

exp(e0/) exp(e0j) exp(-e0/) exp(-e0j) exp(-e^0([X/, Xj])) 



= exp(e0/) exp(e0j) exp(-e0/) exp(-e0j) exp 
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which, thanks to the commutator of left invariant vector fields, depends explicitly on the 
structure of the symmetry group. Note also that the commutator of two vector fields 
results in a factor of for only one power of in the last exponential. Using the previ- 
ous calculations we can easily convince ourselves that this leads to the correct curvature 
components: 

exp(e0/) exp(e0j) exp(-e07) exp(-e0j) exp(-0([eX/, eXj])) 
= (1 + e2e,,,0}0>'= + 0(e3))(l - e'cfj(l>K + 0(e^)) 
= 1 + e\-cfj(j)'k + ei,uMj)Tk + 0(e=^) = 1 + e^F.Vfc + 0{^) . 

We are now in a position to put all ingredients together in order to arrive at our 
regulated classical expression of the Hamiltonian constraint formulated on the regularizing 
manifold S/ F^. Comparing with Eq. (plD, we see that 

nf^ [N] ■= 4{Ki')-^N e^^^ tr(exp(e(/>/) exp{e(j)j) exp(-e0/) exp{-e(j)j) 

UK 

X exp{-e^(P{[Xj,Xj]j){expe(^K,V}) 
has the correct limit (see Eq. ([T^)) 

n^^\N] = lime-=^Hp)[Ar] 

when the regulator is removed. 



4.3 Quantization 



From now on we can again follow the steps of Ref. |T0[ in order to quantize the constraint: 
exponentials of scalars are replaced by (point) holonomies, the volume is quantized to the 
operator of Ref. PD[ and the Poisson bracket is replaced with {ih)~^ times a commutator. 
The result is 



(13) 



UK 



where hj is the holonomy along the J-th regularizing edge, interpreted as multiplication 
operator, and depends on the symmetry group in the following way 



K=l 



What remains to do is to use the decomposition (0) and the quantization of the Eu- 
clidean part of the constraint in order to quantize the Lorentzian constraint. In the reduced 
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formulation the decomposition reads (for a constant lapse function A^) 



n[N] = 2(1 + 1^) [ d^xN{detq)-hr{[Ka,Kb\[E'',E''])-n^^\N] 



= 2{l + i') [ d^xN^o\det{p'K)\-hi{[kj,kj][p',p'])-n^''\N] 




where we defined =: k}u^, K = J d^xK^Ef = Vo/c}pf and used Eq. ([I0|). Regularized 
as above, we obtain 



= 8(l + .2)M-VoiVe"^tr({exp(e0/),ir}{exp(e0j),ir}{exp(e0x),r}) 



which can again be quantized by replacing the exponentials of scalars by (point) holonomies 
and the Poisson brackets by {ih)~^ times a commutator: 



of the integrated extrinsic curvature. 

A non- vanishing cosmological constant A can be included simply by adding the contri- 



Of course, there are factor ordering ambiguities which we ignored in writing the ex- 
pression above. Furthermore, as written down, the operator is not symmetric which would 
have to be achieved by choosing an appropriate factor ordering. In the full theory these 
issues are unsolved and one purpose of the reduced models discussed here can be to gain 
insights by studying them in a simplified, highly symmetric regime. We will, however, not 
enter this discussion because it requires a detailed study of these models, whereas here we 
are mainly interested in general aspects. 

It is immediate to see that the operator ([l^) is manifestly independent of the regulator e 
which was used to reformulate the classical expression in terms of holonomy variables. This 
is so because, in contrast to the classical procedure where the scalar fields are fixed in the 
limit e ^ in which the regulator is removed, in the quantum theory the scalars diverge 
in this point limit in such a way that smeared holonomies are equal to point holonomies 
independent of e. Recall that the quantum configuration space contains distributional 
scalars which are needed to accomplish this result. Therefore, it does not matter whether 



n[N] = 8z(l + L^){dl)-'VoNe''''ti ([hr, K][hj, K][hK, V]^ - H^^'^N] . (14) 



As in the full theory, we use here the quantization 
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we formulate the operator on the auxiliary manifold S/ F, on one of the shrunk manifolds 
S/ or even in the point limit using exclusively point holonomies. In each of these formu- 
lations we have the same state space and the same action of the Hamiltonian constraint 
operator. This feature, which here arose after smearing the scalar fields to regulate classi- 
cal expressions, substitutes the diffeomorphism invariance which in the full theory is used 
to show that the continuum limit removing the regulator there is trivial in the quantum 
theory. Thus, although formulated in a conceptually different way, the quantization of the 
Hamiltonian constraint in the full theory |TD[ and ours for the Bianchi models share the 
same key properties. But of course our discussion was technically simplified by the fact 
that we only have to regard 6-vertices, whereas in the full theory there can be (and have 
to be taken into account) vertices of arbitrary valence. 

It is immediate to see that the constraints are represented in an anomaly free manner: 
the only non-trivial part is to check that diffeomorphisms commute with the Hamiltonian 
constraint which can easily be seen after recalling that diffeomorphisms, i.e. inner automor- 
phisms on the auxiliary manifold, just move the three edges of homogeneous spin networks, 
whereas the Hamiltonian constraint operator fixes them and only changes their labels and 
the contractor. 



5 Hamiltonian Constraints for LRS and 
Isotropic Models 

Similarly to the treatment of volume operators in Ref. we can use the expression ( p^ 
and insert the rotated holonomies /12 = exp(^A3)/ii exp(— to arrive at the constraint 
operator for LRS models, and hi = exp(|A2)/i3 exp(— IA2), /i2 = exp(— |Ai)/i3 exp(|Ai) 
to arrive at the operator for isotropic models. The resulting operator then only contains 
the holonomy operators hi, h^ for LRS models and h^ for isotropic models, and in addition 
operators coming from the exponentials of A/ which manipulate the insertions (see Ref. 

By construction, all these operators are gauge invariant, but the intermediate action 
of holonomies and the volume operator is on gauge non-invariant states. This means that 
the techniques used for gauge invariant states of isotropic models in Ref. have to 



be generalized so as to deal with non-invariant states. In particular, the volume operator 
derived there is not sufficient to calculate the matrix elements of the Hamiltonian constraint 
for isotropic models. 

In contrast to Bianchi models, where the three scalars 0/ and so the three point 
holonomies hj are independent, for LRS and isotropic models there are restrictions for 
the point holonomies (see Eq. (P) leading to the rotated holonomies above. This implies 
that we can apply the following lemma, which has already been used in Ref. [^], in order 
to simplify the constraint operator which contains a product of several holonomy operators. 

Lemma 1 Let g := exp(y4rj) and h := exp{BTj) with A, B E M., i ^ j be matrices in the 
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fundamental representation of SU{2). Then 



gh = hg + h ^g + hg ^ — tT{gh) . 

Let us illustrate this for the flat isotropic model, which is obtained by restricting the 
Bianchi I model to isotropy. We thus start from the expression (|T3|) with all structure 
constants vanishing, and insert the rotated holonomies hi and h2 in terms of h^. Then two 
different holonomies hj and hj, I J, appearing in the operator are always of a form such 
that Lemma |l| applies and we can insert 

hjhj^ = hj^hj + hihj + hj^hj^ - tT{hJ^hj) 

simplifying the sum in Eq. (p!3|) to 

J2 e""^ {tl{[hK, V]) + tl{hj[hK, V]) + tl{hf[hK, V]) - tl{hfhj) tl{hjhf[hK, V])) . 
UK 

The first trace is symmetric in (in fact independent of) / and J and therefore cancels in 
presence of e^"^^ , and the rest can be written as 

^ e^-^^ {tT{{hj - hf)[hK, V]) - ti{hj'hj) ti{hjh-/[hK, V])) . 

UK 

Up to now the explicit expression for the rotated holonomies hi and /12 in terms of 
hs have not been used, but only the fact that Lemma |I] applies because of the special 
character of rotated holonomies. To arrive at the final expression we have to insert the 
explicit formulae for hi, /i2 which leads to a product of /is-holonomies and exponentials of 
Ai, A2. Here we can again apply Lemma |l| in order to collect all factors of h-^ into a single 
power of h^ (except for the one appearing on the right hand side of the volume operator 
when the commutator is written out) and a holonomy independent factor coming from the 
exponentials of A/. Alternatively, we can right from the beginning apply Lemma ^ to the 
definition of hi, h2- 

hi = exp(fA2)/i3exp(-f A2) 

= exp(f A2) (exp(-f A2)/i3 + exp(f A2)/i3 + exp(-f A2)/i3 ^ - tr(exp(-f A2)/i3)) 
= h3 + /ig ^ + exp(7rA2)/i3 - exp(f A2) tr(exp(-f A2)/i3) 

and similarly for /12 such that h^ appears always at the right hand side. After some final 
applications of the lemma we can collect all the factors containing A/ in a single operator 
which manipulates the insertion, whereas the power of h^ changes the spin of isotropic spin 
networks. 

Once the techniques of isotropic spin networks are generalized to gauge non-invariant 
states, the calculation of matrix elements of the Hamiltonian constraint operator in isotropic 
models can be pursued. However, although there is no new input required in addition to 
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the procedure outlined here, the actual calculation can be expected to be quite cumber- 
some and we do not present any further details. These computations are needed when one 
begins to study specific models and are well suited to be done numerically. In LRS models 
the situation is similar: the essential steps have been outlined, but more work had to be 
done for specific models. 

Because for LRS and isotropic models the diffeomorphism constraint vanishes on gauge 
invariant states, the constraint algebra is represented anomaly free in a trivial way. 



6 Discussion 

In this paper we have presented quantizations of Hamiltonian constraints for some cosmo- 
logical models in the framework of loop quantum cosmology. The strategy was to be as 
close to the full theory of loop quantum gravity as possible in order to be able to compare 
with and to draw conclusions for the full theory. We have seen that this is possible to a 
large extent, but with necessary additional input taking care of the symmetry. In partic- 
ular, homogeneous models are formulated in a single point xq and there is no continuum 
limit which in combination with diffeomorphism invariance plays an important role in the 
regularization of the Hamiltonian constraint in the full theory. We have seen that an exten- 
sion of the point xq to an auxiliary manifold, which already appeared in the regularization 
of derivative operators, provides the additional structure. The Hamiltonian constraint is 
regularized on the auxiliary manifold, the extension of which serves as regulator, i.e. the 
classical expression is recovered if the auxiliary manifold is shrunk to the point xq. But the 
quantization of the regulated expression is independent of the extension, and so we have an 
exact analog of the situation in the full theory, where the continuum limit removes the reg- 
ulator classically and the quantization is regulator independent owing to diffeomorphism 
invariance. 

Whereas the investigation of homogeneous models in loop quantum cosmology will be 
more complicated than and very different from the standard treatment of minisuperspace 
quantizations, our constraint operators are, compared with the Hamiltonian constraint 
operator of the full theory, very similar but slightly simpler due to the symmetry. One 
simplification comes from the special nature of vertices being at most 6-valent. Also, the 
action is simpler because the graph of a spin network being acted on is not changed: no 
new edges are created but only the spins of existing ones and the vertex contractor are 
changed. Note that this implies that the large class of special solutions to the Hamiltonian 
constraint found in Ref. has no counterpart in our cosmological models. Furthermore, 



the fact that the constraint algebra is represented anomaly free is realized rather trivially 
for the homogeneous models because the classical algebra is simpler, and so requiring an 
anomaly free representation is less restrictive as compared to the full theory |jl2|, |2^, p6| . 



The calculation of matrix elements of the constraints for Bianchi models can be done 
along the lines of Ref. whereas the computation for LRS or isotropic models would be 



different because techniques for generalized spin networks are needed there. 

Another difference to the full theory is that we do not have to face the issue of extracting 
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a smooth, semiclassical metric from our states. In the full theory, the quantum states 
represent a distributional metric with a discrete structure which, in order to perform a 
classical limit, have to be superposed in some way to semiclassical states. In a homogeneous 
regime, however, the metric is necessarily smooth and each state represents a homogeneous 
(but nevertheless quantum) metric. But we still have to understand how to interpret 
possible solutions of the constraint in a space-time picture or cosmological language, i.e. 
we have to interpret the dynamical constraint as an evolution equation [|17|, llH]. 
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